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ABSTRACT 

An accretion disk around a rotating magnetized star is subjected to magnetic torques which 
induce disk warping and precession. These torques arise generically from interactions between 
the stellar field and the induced surface currents on the disk. Applying these new effects to 
weakly magnetized (B ~ 10 7 -10 9 G) neutron stars in low-mass X-ray binaries, we study the 
global hydrodynamical warping/precession modes of the disk under the combined influences 
of relativistic frame dragging, classical precession due to the oblateness of the neutron star, 
and the magnetic torques. Under quite general conditions, the magnetic warping torque can 
overcome the "Bardeen-Petterson" viscous damping and makes the modes grow. The modes are 
confined to the inner region of the disk, and have frequencies equal to 0.3 — 0.95 (depending 
on the mass accretion rate M) times the sum of the Lense-Thirring frequency, the classical 
precession frequency, and the magnetically driven precession frequency evaluated at the inner 
disk radius r m . As M increases, the mode frequency is reduced relative to the total precession 
frequency at r m since the mode becomes less concentrated around n n due to the increasing 
viscous stress associated with the large M. Because of this, and because the magnetically driven 
precession is retrograde (opposite to the Lense-Thirring precession) and depends strongly on 
M, the mode frequency can have a non-monotonic dependence on the mass accretion rate. This 
may account for several observed features of low-frequency (10-60 Hz) quasi-periodic oscillations 
(LFQPOs) in low-mass X-ray binaries which are otherwise difficult to explain, such as the 
flattening/ turnover in the LFQPO frequency - M correlation or in the LFQPO frequency - kHz 
QPO frequency correlation (e.g., as seen clearly in GX 17+2). 
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1. Introduction 

Accreting neutron stars (NSs) in low-mass X-ray binaries (LMXBs) exhibit rapid, quasi-periodic 
oscillations (QPOs) in their X-ray fluxes. The so-called horizontal-branch oscillations (HBOs), discovered 
in a subclass of LMXBs called Z sources, are low-frequency QPOs (LFQPOs) with centroid frequencies in 
the range of 15-60 Hz, Q-values of order a few, and rms amplitudes <, 10 % (see van der Klis 1995 for a 
review). One interpretation of the HBOs is the magnetospheric beat-frequency model (Alpar & Shaham 
1985; Lamb et al. 1985) in which the HBO frequency is identified with the difference frequency between the 
Keplerian frequency at the magnetospheric boundary and the spin frequency of the NS. However, over the 
last few years, the discovery of kilohertz (300-1300 Hz) QPOs in ~ 20 LMXBs by the Rossi X-Ray Timing 
Explorer (RXTE) has called into question this interpretation of the HBOs (see van der Klis 2000 for a 
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review). In particular, observations indicate that in most Z sources, kHz QPOs (which often come in pairs) 
occur simultaneously with the HBOs. In several atoll sources (which are thought to have weaker magnetic 
fields and smaller accretion rates than the Z sources), similar LFQPOs have also been found at the same 
time when kHz QPOs occur. While the origin of the kHz QPOs is still being debated, it is convenient (as 
in all current theoretical models; see van der Klis 2000 and references therein) to associate one of the kHz 
peaks with the orbital motion in the inner edge of the accretion disk. On the other hand, the discovery of 
nearly coherent oscillations during X-ray bursts in several LMXBs establishes the spin frequency of the NS 
to be ~ 300-600 Hz (see Strohmayer 2000). Obviously, the beating between the kHz orbital frequency and 
the spin frequency cannot produce the 10-60 Hz LFQPOs. 

Stella & Vietri (1998) suggested that the LFQPOs are associated with the Lense-Thirring precession 
of the inner accretion disk around the rotating NS. If the LFQPO and the kHz QPOs are generated at the 
same special radius in the disk, this implies a quadratic relation between the LFQPO frequency and the 
kHz QPO frequency, in rough agreement with observations (usually for certain range of X-ray fluxes; but 
see §4.4). For reasonable values of NS moment of inertia, the data require that the observed LFQPO peaks 
correspond to 2 or 4 times the Lense-Thirring frequency (Stella & Vietri 1998; Ford & van der Klis 1998; 
Psaltis et al. 1999); this may be possible since a precessing disk may induce X-ray flux variation mainly at 
the harmonics of its fundamental precession frequency and there may be "sub-harmonic" feature buried in 
the X-ray power spectra (see Jonker et al. 2000). 

There are several theoretical and observational puzzles associated with the Lense-Thirring interpretation 
of LFQPOs. First, it is well-known that the combined effects of differential precession and viscosity tend to 
keep the inner region of the disk [within ~ (100— 1000)GM/c 2 ] perpendicular to the NS spin axis (Bardeen 
& Petterson 1975). A mechanism to drive the disk tilt is then needed for the precession to make sense. 
Second, because different "rings" in the disk have different precession rates and are strongly coupled to each 
other, a global mode analysis is needed to determine the true precession rate of the inner disk. Markovic & 
Lamb (1998) carried out such an analysis but found that all the modes are damped even when the radiation 
driven warping torque (Pringlc 1996) is included. Third, pure Lense-Thirring precession has difficulty in 
explaining the observed behavior of the LFQPO frequency as a function of the accretion rate of the system 
(Wijnands et al. 1996; Psaltis et al. 1999; see §4.4). 

In this paper, we show that the magnetically driven warping instability (Lai 1999), resulting from 
interactions between the weakly magnetized (<~ 10 8 G) NS and the inner accretion disk (see §2), can 
naturally overcome the Bardeen-Petterson viscous damping and therefore induces a tilt in the inner 
disk. By carrying out a global analysis of disk precession/ warping modes, including magnetic torques 
in addition to the Lense-Thirring and classical precession torques, we obtain, for the first time, growing 
warping/precession modes of the disk (see §3 and §4). Moreover, we show that the magnetically driven 
precession torque (Lai 1999; see §2) significantly affects the actual disk precession frequency as the accretion 
rate increases; this can potentially explain the observed M-dependent behaviors of the LFQPOs (§4). 

2. Magnetically Driven Warping and Precession 

It is known that if the accretion disk is slightly tilted from the equatorial plane of a rotating neutron 
star, it precesses due to general relativistic frame-dragging effect (Lense-Thirring precession) and also 
due to oblateness of the rotating star (classical precession). For a magnetized neutron star, there is an 
additional magnetically driven precession and warping of the disk; these arise from interactions between 
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the stellar magnetic field and the disk before the latter is truncated at the magnetospheric boundary (Lai 
1999). The magnetic torques appear under the generic condition that the stellar magnetic field (treated as 
a dipole in this paper) is not aligned with the stellar spin axis. Depending on how the disk responds to 
the stellar field, two different kinds of torque arise: (i) If the vertical stellar magnetic field B z penetrates 
the disk, it gets twisted by the disk rotation to produce an azimuthal field AB^ — T(B Z that has different 
signs above and below the disk (£ is the azimuthal pitch of the field line and depends on the dissipation 
in the disk), and a radial surface current K r results. The interaction between K r and the stellar B^, gives 
rise to a vertical force. While the mean force (averaging over the azimuthal direction) is zero, the uneven 
distribution of the force induces a net warping torque which tends to misalign the angular momentum of 
the disk with the stellar spin axis, (ii) If the disk does not allow the vertical stellar field (e.g., the rapidly 
varying component of B z due to stellar rotation) to penetrate, an azimuthal screening current will be 
induced on the disk. This interacts with the radial magnetic field B r and produces a vertical force. The 
resulting precessional torque tends to drive the disk into retrograde precession around the stellar spin axis. 

In general, both the magnetic warping torque and the precessional torque are present. For small disk 
tilt angle [3 (the angle between the disk normal and the spin axis), the warping rate and precession angular 
frequency at radius r are given by (see Lai 1999) 

r ^=wlW) cos20 ' (1) 

and ^ 

°™ M = " 2 7rv ?w \n( \ sin2 °' ( 2 ) 
ir z r< \l\r)2^[r)D\r) 

where fi is the stellar magnetic dipole moment, 6 is the angle between the magnetic dipole axis and the spin 

axis, is the orbital angular frequency, and E(r) is the surface density of the disk. The dimensionless 
function D(r) is given by 

D(r) = max (^/r 2 /rf n - 1, ^2H(r)/r in ), (3) 

where H (r) is the half-thickness and r; n is the inner radius of the disk. 

Here we apply these magnetic torques to LMXBs containing weakly magnetized (B ~ 10 8 G) NSs. We 
parametrize the properties of the inner disk using the "inner region" (radiation- and-scattering dominated) 
solution of a-disk (Shakura & Sunyaev 1973)[], with 

E(r) = (1050 g cm-^azlM^M^r^Jir)- 1 , (4) 
H(r) = (1.1 km)M 17 J(r), (5) 

where a = (0.1)a_i is the a-viscosity, M — {1AMq)Mxa is neutron star's mass, M = (10 17 g s _1 )Mi7 is 
the mass accretion rate and r$ — rj (10 6 cm). The dimensionless function J{f) is given by 

(G) 



Such a disk is unstable against thermal- viscous perturbations (see Frank et al. 1992), although the situation is not clear 
when there arc magnetic fields threading the disk. Note that our main results related to the dimensionless mode frequencies 
(see §4.1 and §4.2) are not very sensitive to disk models (see Shirakawa & Lai 2001, where general power-law disks are studied). 

2 Magnetic fields threading the disk can modify J(r) in a model-dependent way (see Lai 1999 for an example). However, the 
basic feature can still be approximated by eq. (ra). 
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where £ is a dimensionless parameter with < £ < 1. Substituting equations (^) and (||) into equations (Q) 
and (0), and assuming the Keplerian flow for f2(r), we get 



r m (r) = (55.3 s- x )C cos 2 6a. lf i 2 26 M 17 J{r) r 6 7 , 
flmW = (-70.4 s- 1 ) sin 2 to.i^M^flfrJ-'jWr- 7 , 

where /i 2 6 = ^/(10 26 G cm 3 ). 

The Lense-Thirring precession angular frequency is given by |^| 

2Gi 3 n s 



(7) 
(8) 



(280 s- L )I 45 



300 Hz 



(9) 



where 1% = (10 45 gcm 2 )/45 is the moment of inertia around the spin axis and Sl s — 2ttv s is the spin angular 
frequency. The classical precession angular frequency is given by 



Q cl (r) 



3G(I 3 ~h) cos/? 



2r 3 



VGMr 



(10) 



where I\ is the moment of inertia around the axis perpendicular to the spin. If we approximate the 
rotating NS as a compressible Maclaurin spheroid (Lai et al. 1993), we have I3 = (l/5)K n M(af + a 2 ,) 
and 1\ = (l/5)K n M(a,2 + a§), where a% is the semi-major axis of the spheroid, and K n is a dimensionless 
parameter which depends on the polytropic index (e.g., «n=l, 0.81482, 0.65345 for n—0, 0.5, 1, respectively). 
For slow rotation, with Cl s = fl s / (GM/ R 3 ) 1 / 2 <C 1 (where R is the radius of a nonrotating star of the same 
mass), the semi-major axis is given by = R(l + q), with (Lai et al. 1994; Appendix A) 



61 = e 2 
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Equations (|J) and ( |l(i| ) then become (with cos f3 ~ 1) 



1 



l+2«r. 



1 

5 



n iT (r) = (313 s" 1 )/^ 
ft*) - (- 1 95 S ->K(1-^)( 55 ^_)'M,- 1 "W 



300 Hz 

-7/2 



6 ' 



(11) 

(12) 

(13) 
(14) 



where i?6 = i?/(10 6 cm). Morsink & Stella (1999) have given exact numerical solutions of Qlt and f2 c i 
for a variety of nuclear equations of state. But analytic expressions ( p^ ) and (14) provide the dominant 
contributions, and are adequate for the purpose of this paper. 



3. Global Precession/Warping Modes: Equations 

Since the precession rates FIlt, ^ci, depend strongly on r, coupling between different rings is 
needed to produce a global coherent precession. Such coupling can be achieved through either viscous 



3 In the strong gravitational field and/or rapid stellar rotation regime, the Lense-Thirring precession (angular) frequency is 
replaced by the nordal precession frequency Q noc i = Q — Qg with Q the orbital frequency and Qg the vertical frequency. The 
fractional deviation of fi no d from Qlt is of order ^(M/r) 1 / 2 (where a = 70 3 /M 2 ) and will be neglected in this paper. 
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stress or bending waves (e.g., Papaloizou & Pringle 1983; Terquem 1998 and references therein). In the 
viscosity dominated regime, the formalism of Papaloizou & Pringle (1983) can be used (see also Pringle 
1992; Ogilvie 2000; Ogilvie & Dubus 2001; note that in the linear regime, the different formalisms of disk 
warping dynamics are equivalent). We specify a warped precessing disk by the disk normal vector l(r, t). 
In the Cartesian coordinate, with the z-axis along the NS spin, we write 1 = (sin (3 cos 7, sin (3 sin 7, cos j3), 
with /3(r,t) the tilt angle and 7(7", t) the twist angle. For (3 <C 1, the dynamical warp equation for 1 (see Lai 
1999) reduces to an equation for W(r, t) = /3(r, £)e ry ( r ' t ): 



dW 



3/y 2 
At 



1 



It J' 



2r 



dW 
dr 



1 d 2 W 



+ i {Sl LT + n d + Sim) w + T m W, (15) 



where is the viscosity which tends to reduce the disk tilt. We assume that the ratio of to the usual 
viscosity v\ is constant. For a disk, v\ = aH 2 Sl, the viscosity rate is 



3 (r) 



v 2 {r) 



(16.7 s" 1 ) 



(16) 



To look for global modes we consider solutions of the form W(r,t) = e l(7t W(r) with the complex mode 
frequency a (= a r + ie>i). It is convenient to define a dimensionless mode frequency a = <7T v i sc (ri n ), where 
TVisc^in) is given by eq. ( |l6| ) evaluated at r- m , the inner disk radius. Similarly, we define dimensionless 
quantities: 



SIlt 

Slcl 
f 

J- m. 



16.8 145 ( Vs ) 
5 V300 Hz/ 

11.7 k 2 , ( 1 



f^LT (nn)Tvisc(nn) = 
^cl(^*in)^~visc (^*in ) ~~ 

^m(n n )T visc (r in ) = -4.2 sin 2 ^ 
r m (r in )r v i sc (r in ) = 3.3Ccos 2 6» 



- 1 M~ 1/2 M~ 2 7- 2 ( T[n \ 
-1^1.4 M 17 J. m ^ iQkm j 



1/2 



n 
" 5 

V-2 



a-_\M^lM^Rljr 2 , 



f-VV 

V300 Hz/ 

7/2 



-7/2 



(17) 
(18) 
(19) 
(20) 



where D in = D{r- ln ) and J ln = J(r\ n ). Note that we can use eq. ( |13[ ) to obtain more explicit expression of 
SIlt- Equation ( |l5| ) can now be reduced to the dimensionless form: 



iaW 



4a .5/2 

_ 1 J 2 d 2 W 
~ 2x 3 / 2 J?~ck? 

where x = r/r- 1VL and J' = dj jdx. 
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(21) 



It is clear from eq. ( |2l| ) that <r depends only on five dimensionless parameters SIlt, Sl c \, Sl m , T m , and 
v\jvi as well as two dimensionless functions D{x) and J{x). To obtain a in physical units we need to know 
7*i n . We adopt the simple ansatz: 

r m = max(r m , ri SC o ) , (22) 



where the magnetosphere radius r m is given by 



18?7/J26 7 M 14 1/7 M 17 2/7 km, 



(23) 
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(with r\ ~ 0.5), and the inner- most stable circular orbit Q nsco is given by 

r ISCO = 6GM/c 2 = 12.4 M 1A km. (24) 

The critical mass accretion rate Myj^ c below which r m = r m and above which r- m = nsco is obtained by 
equating eqs. (^) and ( |j) as: 

Mi7, c = 3.7 r^^M^. (25) 

A more elaborate prescription for the inner disk radius when both general relativity and the magnetic field 
are important is discussed in Lai (1998). 

To solve eq. (|2l]) for the complex eigenfunction W(x) and eigenvalue a, six real boundary conditions 
are needed. In our calculation, the disk extends from x- m = 1 to x ou t = 50. For large x and large |a|, 
equation ([H]) can be solved analytically, giving 



W(x) ex exp 



(26) 



where we should choose the sign of (itr) 1 / 2 so that W(x) — > as x — > oo. This approximate analytical 
solution, evaluated at x QU t, together with its derivative, gives four (real) outer boundary conditions. The 
inner boundary condition generally takes the form W'(x- m ) — aW(xi n ), with a being a constant. Most 
of our results in §4 will be based on a = (corresponding to zero torque at the inner edge of the disk), 
although we have experimented with different a's (see Figs. 1-2) and found that for \a\ <, 1 our results are 
unchanged to the extent that there are other parameters in the problem which have greater uncertainties. 
In numerically searching a mode, we make a guess for the eigenvalue a and integrate eq. (|l]) from x out to 
X[ n using the Kaps-Rentrop scheme (Press et al. 1992). We use the globally convergent Newton method to 
find the correct value of a that satisfies the boundary conditions. 



4. Numerical Results and Discussion 

In this section, we first study numerical properties of equation (||l]) and then discuss specific cases 
relevant to accreting NSs in LMXBs. 



4.1. Mode Eigenfunction and Eigenvalue 

For a given set of parameters (&lt, f^ci, fj mi T m ), equation ( |2l| ) allows for many eigenmodes. Here we 
shall focus on the "fundamental" mode which is more concentrated near the inner edge of the disk and has 
larger oy (global precession frequency) and smaller &i (damping rate) than any other "higher-order" modes 
(see Fig. 1). 

If a mode were infinitely concentrated at the inner radius of the disk, one expects that the mode 
frequency ay is just the sum of the frequencies evaluated at the inner disk radius, i.e., ay = Q,lt + ^ci + ^m- 
However, the calculated ay is always smaller than CIlt + ^ci + because the mode is not infinitely 
concentrated but has a finite width. 



4 The correction to nsco due to stellar rotation is negligible since a ~ 0.1 for NSs in LMXBs. 
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Figure 1 shows the tilt angle 0(x,t = 0) — \W(x)\ associated with the modes for different sets of 
(OiTjfm). We have set £l c \ = f2 TO = (since they play a similar role as CIlt), J{%) = 1, and v\/v2 = 1 in 
eq. for simplicity^. We see that as (Ilt an d f m increase, the fundamental modes (solid lines) become 
more concentrated near the inner radius of the disk. This behavior can be understood heuristically: for a 
given £Ilt(t), a larger VLlt implies smaller viscosity, and thus the coupling between different disk radii is 
reduced. 

Figure 2 shows the mode frequency a in units of VI lt ( r in), or a/fi^T^m) = o'/Olt, as a function of 
Qlt for different values of T m [We again set Vl c \ = O m = 0, J(x) = 1, and vi/i>2 = 1]. We see that & r /ClLT 
always lies between 0.4 to 0.95 for the relevant ranges of n LT (10 to 10 4 ). The ratio a r /ft LT increases and 
approaches unity as Qlt and T m increase. This is consistent with the behavior of the mode eigenfunction 
(see Fig. 1) that a larger VLlt or T m makes the mode more concentrated near the inner disk edge. 

4.2. Global Warping Instability Criterion 

As discussed in §1, differential precession tends to be damped by the v-i viscosity. Figure 2(b) shows 
that di/VtLT — Pi/^LT (f*in) is always positive (which implies damping) for T rn = and lies between 0.25 
(for CIlt = 10) and 0.04 (for CIlt = 10 ), corresponding to Q value (Q w ay/o^) between 2 and 20 (the 
high-Q regime was also explored by Markovic & Lamb 1999 where a different disk model was adopted). 
We see that ay decreases as T m increases, and becomes negative (implying mode growth) when the ratio 
r m /riiT is sufficiently large. The numerical values for &i can be approximated by 

a t =-af m + bn° L 7 T , (27) 

with a ~ (0.5 — 1.0) and b ~ 0.5. For the mode to grow (ay < 0) we require 

f m £ Q°t Global Warping Instability. (28) 

For Mi 7 < Mi7 ;C [see eq. (^5|)1, so that rj n = r m , this condition becomes 

0.7< rt« eatll^MUMif (^) "' 3 ( JL) M * 1. (29) 

For Mi 7 > Mi7 iC , so that r in = riscOi we require 

O^CMleCos^V-^MiV (_^_)"°- 7 ^y' 3 £ l. (30) 

We see that for parameters that characterize accreting NSs in LMXBs the mode growth condition can be 
satisfied (see §4.3 for specific examples), although not always. In general, high (but not unreasonable) ( (> 
a few) and M17 are preferred to obtain growing modes. 

4.3. Dependence of Mode Frequency on M 

In §4.1 and §4.2, we have seen how SIlt and r m affect the dimensionless mode frequency a (— ay + 
while setting O c i = f2 m = 0. We have found that ay is smaller than the naively expected value &lt for small 



5 Note the results shown in Figs. 1-2 are not sensitive to J since J-^ has been absorbed in the definitions of dimensionless 
frequencies, eqs. (^t|)-(|2o[). See also §4.3. 
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CIlt and f m , but becomes closer to CIlt as O^r and f m increase. We have also found that a sufficiently 
large T rn tends to make the mode unstable. In this subsection, we consider some specific examples to 
illustrate the dependence of the mode frequency a on M and other parameters for NSs in LMXBs. 

Figure 3 shows the global precession frequency (i.e., a r /2ir] the real part of the mode frequency in 
physical units) Q as a function of the mass accretion rate for a fixed parameter set [Parameter Set A: 
Mi >4 = l,Re = l,n = l,v a = 300 Hz, a_i = 1, C = 5, sin 2 9 = 0.1, v 2 /v 1 = 1, = 2, r? = 0.5, £ = 0, 
and D(x) as given in eq. (§)], for three different cases: (i) only the Lense-Thirring precession is included, 
(ii) both the Lense-Thirring and classical precessions are included, and (iii) the magnetic precession and 
warping are included in addition to the Lense-Thirring and classical precessions. It also shows the fiducial 
precession frequency, z^d, which is the sum of the considered frequencies evaluated at the inner radius of 
the disk, i.e., v Rd = (O ir (r in ) + ft c i(r in ) + £l m (r in ))/27r. 

The dependence of z/fid on M is understood in the following manner. From eqs. (|T3|), (14), and (J2J) , 
together with eqs. (|2^)-(|24|), we find fij/r^in) oc A/ 6 / 7 , f2 c i(rjn) oc M, and fl m { r in) oc M 3 when r- m — r m , 
and f^x^in) oc M , f2 c l(?"in) oc M°, and f2 m (ri n ) oc M when r; n = nsco Thus we obtain, for the above 
three different cases: 

(i) iv fid cx G(M 6 / 7 ) (Mi 7 < Mi 7)C ), VRA cx 0(M°) (M i7 > M ir , c ). (31) 

(ii) i/fid cx 0(M 6 / 7 ) - O(M) (Mir < Mi7,c), ^fid cx 0(M°) (M 17 > M 17 , c ). (32) 

(iii) z/ fid oc G(M 6 / 7 ) - O(M) - 0(M 3 ) (Mi 7 < % c ), (33) 
z/ fid oc 0(M°) - O(M) (Mi 7 > M i7 , c ). (34) 

For Parameter Set A adopted in Fig. 3, the transition from r ul — r rn to r- m — nsco occurs at Mi7 iC = 1.3 
[see eq. (p5|)]. 



We see from Fig. 3 that the behavior of the mode frequency a r /2w as a function of M is different from 
that of z/fid for all three cases (i)-(iii). For case (i) and (ii), a r /2TT is as small as ~ 30% of i/fid for high M 
(M17 ~ 1.5 — 2.0) and the dependence of a r /2ir on M is very different from that of v&&. This is understood 
from our discussion in §4.1: The viscosity rate increases (and &lt + H c i decreases) as M increases [see 
eq. (O)]; this tends to spread the mode by coupling (through viscous stress) disk rings at different radii, 
thereby reducing a r /2ir relative to z^d- The increasing importance of viscosity for high mass accretion rate 
is also expected for case (iii), where additional magnetic effects are included. However, we see from Fig. 3 
that a r /2n becomes closer to fgd when the magnetic warping torque (r m ) is included. This is understood 
from our finding in §4.1 (see Figs. 1-2) that the presence of r m makes the mode more concentrated near 
the disk inner radius. 

Figure 4 shows the mode frequency a r /2ir as a function of M for three sets of parameters, Set A, B, 
and C, where ((U26; sin 2 9) = (2, 0.1), (2, 0.5) and (4, 0.1), respectively, while the other parameters are fixed 
to the standard values: M1.4 = 1, Re = 1, n = 1, v s = 300 Hz, a_i = 1, £ = 5, 1^2/^1 = 1> f] = 0.5, £ = 0. 
All of the Lense-Thirring, classical precessions and magnetic precession and warping are included [case 
(hi)]. For Parameter Set A and B, the inner radius of the disk reaches nsco at the same Mn^ c = 1.3 [see 
eq. (p5|)]. For Parameter Set A, |ri m (r in )| is much smaller than f2z,T(?~m) + ^ci^in) for small M17 although 
it becomes more important for large M17; thus at M17 = Mi7 jC , <r r /27r suddenly changes from an increasing 
function of M as given by eq. (^3|) (with small second and third terms) to a linearly decreasing function of 



6 In Figs. 3-6, we multiply the mode frequency by 4 to facilitate comparison with observations, since the QPOs appear to 
manifest as harmonics of the fundamental frequency (Stella & Vietri 1998; Psaltis et al. 1999; see also §5). 
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M as given by eq. (0). For Parameter Set B, |fl m (7m)| K sm2 ^ is 5 times larger than that in Parameter 
Set A and is important even for small Mvj < Mn^ c ; thus a r /2n shows a gradual turnover feature according 
to eq. (HH) and becomes a linearly decreasing function of M [see eq. above M 17 c . For Parameter Set 
C, the inner radius of the disk is always given by r m for the relevant M (since M17 c = 5.2). As |f2 c i(?"in)| 
and |0 TO (ri n )| are much smaller than fiz,T(nn), a r /27r is a monotonically increasing function as given by 
eq. ( |33| ) (with small second and third terms). 

Figure 5 illustrates the dependence of the mode frequency a r /2n on the polytropic index n and the 
dimensionless parameter £ [defined in eq. (^|)] (other parameters are the same as in Parameter Set A 
adopted in Fig. 4). Most nuclear equations of state have effective polytropic index between 0.5 and 1. We 
see that a r /2ir is larger for a stiffer equation of state (smaller n) and a larger £, although the dependences 
are not as strong as on other parameters (cf, Fig. 4). 

4.4. Comparison with the Observed Behaviors of LFQPOs 

The behavior of oy as a function of M (see Figs. 3-5) is similar to the features observed for the 
LFQPOs in LMXBs (see van der Klis 1995, 2000). While for some systems, the LFQPO frequencies increase 
as the X-ray flux increases, for other systems, a non-monotonic (even opposite) correlation between the 
QPO frequency and M has been observed. As an example, for the Z source GX 17+2, using the position 
of the source on the X-ray color-color diagram as an indicator of M, Wijnands et al. (1996) found that 
the HBO frequency increases with M for small M but decreases with further increase of M . In the simple 
model where the LFQPO is identified with the precession frequency at the disk inner radius r- m (Stella & 
Vietri 1998; see §1), the non-monotonic behavior is difficult to explain, even when the classical (retrograde) 
precession is included^] (see Psaltis et al. 1999; Morsink & Stella 1999). However, we see from Fig. 4 that 
when we consider global disk modes and include the magnetic effects, a variety of ay — M behaviors are 
possible. In our model the turnover of the oy — M correlation occurs for two reasons: (1) the magnetically 
driven retrograde precession Q m cx — 1/r 7 becomes more important with increasing M [see eqs. ( j33[) & 
(p4|)]; (2) the viscous stress becomes more important as M increases, and makes the mode less concentrated 
around r; n ; thus ay is more reduced from with increasing M. We note that even when magnetic effects 
arc not included [case (i) & (ii)], the "turnover" feature is obtained due to mechanism (2) as explained 
above (see Fig. 3). However, these modes are highly damped by viscosity. Only the magnetic warping 
torque can make the mode grow (see Fig. 4). 

A related puzzle concerns the correlation between the LFQPO frequency and the higher kilohertz QPO 
frequency (e.g., Stella & Vietri 1998; Ford & van der Klis 1998; Psaltis et al. 1999; see §1). It was found 
that for several sources the positive correlation flattens (and even turns over) for large kHz QPO frequencies 
(and large M) (see Homan et al. 2001 for a clear case of such turnover in GX 17+2). Again, this feature 
is difficult to explain with the simple picture of Lense-Thirring and classical precessions at n n , but can be 
accounted for qualitatively by our global disk oscillation model with magnetic effects (see Fig. 6; note that 
we have not tried to vary the model parameters to "fit" the observations). More quantitative comparison 
between our models and the data require a better understanding of the origin of kHz QPOs (e.g., it is not 
certain whether the kHz QPO frequency is the Keplerian frequency at r; n ; see §5 for other uncertainties 
expected in real systems). 



7 This is because Qlt oc r 3 and O c i oc — r 7 / 2 have a similar r-dependence; see eqs. (113) and Q14I 
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5. Concluding Remarks 

We have shown in this paper that the inner region of the disk around a weakly magnetized (~ 10 8 G) 
NS is warped due to magnetic field - disk interactions (see Lai 1999) and therefore must precess under 
the combined effects of relativistic frame dragging, classical precession (due to the oblateness of the NS) 
and magnetic torques. We have found growing warping/precession modes of the inner disk and shown that 
these modes have properties that resemble the 10-60 Hz low- frequency QPOs observed in LMXBs (see §4). 

Although our treatments of disk warping/precession go beyond many previous models of QPOs 
(these models typically identify a QPO frequency as the characteristic frequency of a test mass at certain 
special radius in the disk; see van der Klis 2000 for a review), they still contain important idealizations. 
For example, we have assumed the stellar field to be dipolar when considering the magnetic field - disk 
interactions. This is unlikely to be correct for LMXBs since the disk (with r; n only a few stellar radii) can 
produce significant change in the global field topology even when the intrinsic stellar field is dipolar (see Lai, 
Lovelace & Wasserman 1999). Observationally, the absence of persistent millisecond pulsation in all known 
LMXBs but SAX J1808. 4-3658 (Wijnands & van der Klis 1998) and the frequency drift of burst oscillations 
(e.g., Strohmayer 2000; see also Cumming et al. 2001) may imply that the NS does not possess a well-defined 
dipole field, although a more complex field is possible. Also, we have only studied the linear behavior of the 
warping/precession modes in this paper. Because of these idealizations and the uncertainties in observations 
(e.g., it is difficult to determine M precisely from observations), more quantitative comparison between the 
observational data and our theory is premature at present. Nevertheless, the results presented in this paper 
demonstrate that magnetically driven warping and precession can give rise to a variety of new possibilities 
for the dynamical behaviors of inner accretion disks around magnetic NSs. 

Concerning the observed properties of QPOs, many issues remain unanswered in our paper (and in 
other related QPO studies): e.g, How does the mode manifest itself as a variation in the X-ray flux? How is 
the observed QPO amplitude (as large as 15%) produced? One possibility is that the LFQPO is caused by 
occultation of the central NS by the warped, precessing inner disk. Another issue is that for the observed 
NS spin rates (based on observations of burst oscillations) and reasonable moment of inertia of the neutron 
star, the data require that the observed LFQPO corresponds to 2 or 4 times the precession mode frequency 
(see §1; the fact that the global mode frequency is smaller than the precession frequency at r- m exacerbates 
the problem); while this may be possible, more theoretical study on this issue is needed. Finally, there is 
an indication (but by no means proven) that some black hole systems display low-frequency QPOs which 
resemble their counterparts in NS systems (Psaltis, Belloni & van der Klis 1999). Clearly, the magnetic 
warping and precession torques discussed in this paper are irrelevant to black hole systems^]. We note 
that the phenomenology of black hole QPOs is much less developed, and that many of the Ai-dependent 
behaviors of LFQPOs observed in NS systems are absent in the black hole systems (see Remillard 2001). 

This work is supported in part by NSF Grant AST 9986740 and NASA grant NAG 5-8484, as well as 
by a research fellowship (to D.L.) from the Alfred P. Sloan foundation. 



8 Note that even without the magnetic torques and the classical precession, the mode frequency is reduced relative to the 
Lense-Thirring frequency at r; n as M increases (see Fig. 3); this result applies equally to the black hole systems except that 
the mode is damped in the absence of other excitation mechanisms. 
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Fig. 1. — The disk tilt angle of the warping/precession modes for various parameter sets. The solid curves 
represent the fundamental modes for (£l LT , f m ) =(1000, 0), (100, 100), and (100, 0) (with Cl c i = Q m = 0) 
from left to right, with the corresponding mode frequency a = (<3y,<7i) = (860,78), (80,-46), and (70,15); 
the dotted curve represents a higher order mode for (CIlt, r m )=(100, 100), with a = (a r ,&i) = (41,6.0); all 
these are calculated using the inner boundary condition (B.C.) W = 0. The dashed curve shows the case for 
(CIlt, T m ) =(100, 100), with the mode frequency a = (65, —37), calculated using the inner B.C. W = W. 
The eigenfunctions are normalized such that the maximum tilt angle is 1. 
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Fig. 2. — The upper panel shows the mode frequency ay in units of Qlt as a function of Qlt for different 
values of r m (with O c i = Q m = 0). The lower panel shows the corresponding mode damping rate <7i (in 
units of VIlt)- Note that negative &i implies growing mode. The dot-dashed curves are obtained using the 
inner boundary condition W — W, while all other curves are based on B.C. W = 0. 
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Fig. 3.— The global precession mode frequency a r /2iT as a function of the accretion rate Myj. The 
frequencies are multiplied by 4 to facilitate comparison with observations (see the text for discussion). 
The upper three curves represent the fiducial frequency z/gd for three different cases described in the text: (i) 
only the Lense-Thirring precession is included (long dashed curve), (ii) both the Lense-Thirring and classical 
precessions are included (short dashed curve), and (iii) the magnetic precession and warping are included as 
well as the Lense-Thirring and classical precessions (solid curve). The thick solid curve is the corresponding 
global precession frequency a r /2n for case (iii) (Note that the modes grow above Myj ~ 0.6; marked by a 
dot in the figure). The lower two dashed curves are the corresponding global precession frequency oy/27r for 
case (i) and (ii). The lower light solid curve shows a r /2n for case (iii) obtained with the magnetic warping 
torque artificially turned off. 
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Fig. 4. — Correlation between the mode frequency a r /2w and the accretion rate M\-j for Parameter Sets A, 
B, and C: (^26, sin 2 9) = (2, 0.1), (2, 0.5), & (4, 0.1), respectively (see the text for the values of other fixed 
parameters) when all of the Lensc Thirring, classical, and magnetic precessions and magnetic warping are 
taken into account. The solid portion of the curve corresponds to growing mode and the dashed portion 
corresponds to damping mode. To compare with LFQPOs observed in LMXBs, we multiply a r /2-K by 4. 
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Fig. 5. — The mode frequency a r /2ir (multiplied by 4) as a function of Mn for polytrope index n = 0.5 
(upper two curves) and 1 (lower two curves), and for £ = (solid curves) and £ = 0.5 (dashed curves) [see 
ec l- ©]■ Other parameters are the same as in Parameter Set A adopted in Fig. 4. 
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Fig. 6. — Correlation between the mode frequency a r /2ir (multiplied by 4) and the Keplerian frequency at 
the inner radius of the disk for Parameter Sets A, B, and C as in Fig. 4. The solid portion of the curve 
corresponds to growing mode and the dashed portion corresponds to damping mode. 



